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ABSTRACT 

XTE J1946+274 = GRO J1944+26 is a 15.8 s Be/X-ray pulsar discovered simultaneously 
in 1998 September with the Burst and Transient Source Experiment (BATSE) on the Compton 
Gamma Ray Observatory (CGRO) and the All-Sky Monitor (ASM) on the Rossi X-ray Timing 
Explorer (RXTE). Here we present new results from BATSE and RXTE including a pulse timing 
analysis, spectral analysis, and evidence for an accretion disk. Our pulse timing analysis yielded 
an orbital period of 169.2 days, a moderate eccentricity of 0.33, and implied a mass function of 
9.7 Mq. We observed evidence for an accretion disk, a correlation between measured spin-up 
rate and flux, which was fitted to obtain a distance estimate of 9.5 ± 2.9 kpc. XTE J1946-I-274 
remained active from 1998 September - 2001 July, undergoing 13 outbursts that were not locked 
in orbital phase. Comparing RXTE PGA observations from the initial bright outburst in 1998 
and the last pair of outbursts in 2001, we found energy and intensity dependent pulse profile 
variations in both outbursts and hardening spectra with increasing intensity during the fainter 
2001 outbursts. In 2001 July, optical Ha observations indicate a density perturbation appeared 
in the Be disk as the X-ray outbursts ceased. We propose that the equatorial plane of the Be 
star is inclined with respect to the orbital plane in this system and that this inclination may be 
a factor in the unusual outburst behavior of the system. 

Subject headings: accretion—stars:pulsars:individual:(XTE J1946-I-274, GRO J1944-I-26)—X-rays: 
binaries 
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1. Introduction 

Be/X-ray binaries are the most common type of 
accreting X-ray pulsar systems. They consist of a 
pulsar and a Be (or Oe) star, a main sequence star 
of spectral type B (or O) that shows Balmer emis¬ 
sion lines (See e.g., Slettebak 1988 and Coe 2000 
for reviews.) The line emission is believed to be 
associated with an equatorial outflow of material 
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expelled from the rapidly rotating Be star which 
probably forms a quasi-Keplerian disk near the Be 
star (Hanuschik 1996; Quirrenbach et al. 1997). 
X-ray outbursts are produced when the pulsar in¬ 
teracts with this disk. Be/X-ray binaries typically 
show two types of outburst behavior: (a) giant 
outbursts (or type II), characterized by high lu¬ 
minosities (Lx ^ 10^^ ergs s“^) and high spin-up 
rates (i.e., a significant increase in pulse frequency) 
and (b) normal outbursts (or type I), character¬ 
ized by lower luminosities (Lx ~ 10^® — 10®^ ergs 
s“^), low spin-up rates (if any), and recurrence at 
the orbital period (Stella, White, & Rosner 1986; 
Bildsten et al. 1997). As a population Be/X-ray 
binaries show a correlation between their spin and 
orbital periods (Corbet 1986; Waters & van Kerk- 
wijk 1989). 

On 1998 September 8 a 15.8 s pulsar, which was 
designated GRO J1944-I-26, was discovered with 
the Burst and Transient Source Experiment (Fish¬ 
man et al. 1989, BATSE) on the Compton Gamma 
Ray Observatory (CGRO). GRO J1944-I-26 was 
localized to a 5° x 8° error box (Wilson et al. 
1998). At the same time, the All-Sky Moni¬ 
tor (Levine et al. 1996, ASM) on the Rossi X- 
ray Timing Explorer (RXTE) discovered a new 
source which they localized to a 6' x 26' error 
box and designated XTE J1946-I-274 (Smith & 
Takeshima 1998). Subsequent observations with 
the Proportional Gounter Array (Jahoda et al. 
1996, PGA) also on the RXTE on 1998 Septem¬ 
ber 16 revealed 15.8 s pulsations, confirming that 
BATSE and RXTE were seeing the same object 
(Smith & Takeshima 1998). Scanning observa¬ 
tions with the PGA further improved the position 
to a 2.4' error circle (Takeshima & Chakrabarty 
1998). Observations with BeppoSAX further im¬ 
proved the X-ray position to R.A. = 19M5to38s, 
Deck = -|-27°21'.5 (equinox 2000.0; error radius 
1' at 95% confidence) (Campana, Israel, & Stella 
1998). XTE J1946-I-274 lies in the error box of the 
Ariel 5 transient 3A 1942-1-274; however, there is 
a non-negligible probability that this is a chance 
association (Gampana, Israel, & Stella 1999). 

A likely optical counterpart has been identified 
as an optically faint B ~ 18.6 mag, bright infrared 
(i7 ~ 12.1) Be star at R.A. = 19Ii45m39s, Deck 
= 27°22'00", which lies in the refined BeppoSAX 
error circle(Verrecchia et al. 2002). Spectroscopic 
and photometric data indicate a BO-lIV-Ve star at 


8-10 kpc. The counterpart’s mass is not well con¬ 
strained, but is expected to be in the 10 — IOMq 
range. Two previously proposed candidates (Is¬ 
rael, Polcaro, & Covino 1998; Ghavamian & Gar¬ 
cia 1998) lie 10" outside the BeppoSAX error cir¬ 
cle. 

Using data from the 1998 outburst taken with 
the PGA and the High Energy X-ray Timing Ex¬ 
periment (Rothschild et al. 1998, HEXTE) on 
RXTE, Heindl et al. (2002) revealed a cyclotron 
resonance scattering feature or cyclotron line near 
35 keV. This feature implies a magnetic field 
strength of 3.1(1 + z) x 10^^ G, where z is the 
gravitational redshift of the emission region. Us¬ 
ing about 1 year of RXTE ASM 2-10 keV flux 
measurements, beginning with the 1998 outburst, 
Gampana, Israel, & Stella (1999) observed 5 out¬ 
bursts in which they found evidence of an ~ 80 
day modulation in the outburst flux. They found 
that the first outburst in September 1998 was 
markedly different in rise and decay timescales 
than the 4 following outbursts and based on this 
suggested that the first outburst might be classi¬ 
fied as a giant outburst. XTE J1946-I-274 was ob¬ 
served with the Indian X-ray Astronomy Exper¬ 
iment (Paul et al. 2001, IXAE) 1999 September 
18-30 (MJD 51408-51421) and 2000 June 28-July 
7 (MJD 51723-51733). They detected 15.8 s pul¬ 
sations in the 2-6 and 6-18 keV bands, with similar 
double peaked profiles in both observations. The 
pulse period history was consistent with a constant 
intrinsic spin-up plus an eccentric orbit. The re¬ 
ported pulse periods were P = 15.78801(04) s and 
on MJD 51445.0 and P = 15.76796(18) s on MJD 
51727.5. 

XTE J1946-I-274 was active from 1998 Septem¬ 
ber through 2001 July, when it dropped below 
the detection threshold of the RXTE PGA. As 
of 2002 July, no additional outbursts have been 
seen with the RXTE ASM. Table 1 lists the dates 
and instruments observing XTE J1946-|-274’s lo¬ 
cation. In this paper we will present results from 
observations of XTE J1946-I-274 with BATSE and 
RXTE, including histories of pulse frequency and 
total flux for XTE J1946-I-274. From the pulse fre¬ 
quency history, we derive an orbital solution. Ap¬ 
plying this orbital solution, we will investigate the 
observed correlation between spin-up torques and 
observed flux and discuss implications for the ac¬ 
cretion mechanism. We will describe the unusual 
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orbital phasing of the outbursts. Using data from 
the RXTE PCA and HEXTE, we will investigate 
spectral variations in the initial bright outburst 
and the last 2 outbursts. We compare our X-ray 
results to optical Ha observations. 

2. Analysis and Results 
2.1. Frequency Search 

To determine pulse frequencies for XTE J1946-I-274, 
we performed a grid search over a range of candi¬ 
date frequencies using data from the Large Area 
Detectors on BATSE(Fishman et al. 1989). This 
search technique is described in detail in Finger 
et al. (1999) and Wilson-Hodge (1999). We will 
briefly summarize the technique here. Variations 
of this technique, including searches in pulse pe¬ 
riod, have been widely used (Buccheri et al. 1983; 
Bildsten et al. 1997, for example). 

The pulsar search technique consists of 3 steps 
(1) data selection and combination, (2) 20-50 keV 
pulse profile estimates, and (3) a grid search in fre¬ 
quency. First, the BATSE DISCLA channel 1 (20- 
50 keV, 1 s time resolution) data were selected for 
which the source was visible, the high voltage was 
on, the spacecraft was outside the South Atlantic 
Anomaly, and no electron precipitation events or 
other anomalies were flagged by the BATSE mis¬ 
sion operations team. The count rates were com¬ 
bined over the 4 LADs viewing XTE J1946-I-274, 
using weights optimized for an exponential energy 
spectrum, f{E) = Aexp{—E/kT) with temper¬ 
ature kT = 12 keV, and grouped into « 300 s 
segments. A segment length of 300 s was used 
because it was short enough that the background 
was well-fitted by a quadratic model. The segment 
boundaries were chosen to avoid Earth occultation 
steps from bright sources(Wilson-Hodge 1999). 

The second step in this process was estima¬ 
tion of an initial 20-50 keV pulse profile for each 
segment of data. In each segment, the com¬ 
bined rates were fitted with a model consisting 
of a sixth-order Fourier expansion in pulse phase 
(representing the 20-50 keV pulse profile) r{t) = 
a/ifc exp(t27rIi())o(t)) where ahk is the esti¬ 
mated complex Fourier coefficient for harmonic^ h 
in segment k and a spline function with quadratics 

^In this paper, harmonics are defined as nu where n = 

1, 2, 3,... and u is the pulse frequency. 


in time (representing the background plus mean 
source count rate). A sixth-order Fourier expan¬ 
sion was chosen based on the number of har¬ 
monics that are significant in a 1-day observa¬ 
tion, the typical integration time required to de¬ 
tect XTE J1946-I-274. Our initial pulse phase 
model was of the form (fioit) = h'o{t — to), where 
t'o = 63.1513 mHz before MJD 51270 and i^o = 
63.2800 mHz after MJD 51270 were barycentric 
frequencies and t was the time corrected to the 
Solar System barycenter® using the JPL DE-200 
ephemeris (Standish et al. 1992). The initial pulse 
frequency was changed after MJD 51270 because 
continued spin-up of the pulsar moved the ob¬ 
served frequency too far away from the folding fre¬ 
quency. The value and slope of the spline function 
were required to be continuous across adjacent 300 
s segment boundaries, but not across data gaps. 

Ideally we would like to be able to fit a pulse 
profile to the entire 4-day interval of data for each 
frequency grid point; however this is computation¬ 
ally very very expensive. Instead we have formu¬ 
lated an equivalent method in which we first fit 
pulse profiles to short segments of data using a 
single fixed frequency. Then the profiles are com¬ 
bined and shifted by a frequency offset correspond¬ 
ing to each grid point. Fitting pulse profiles to 
short data segments has other advantages in ad¬ 
dition to improving computational efficiency. (1) 
If the folding frequency chosen was incorrect, the 
profile will be smeared out. However, the degree 
of smearing in pulse phase A(j) is approximately 
given by A(j) = AvEgn, where Av is an offset in 
frequency. Eg is the segment length, and n is the 
number of Fourier coefficients. For Eg = 300 s, 
Av = 0.17 mHz, and n = 6, A(f> = 0.31, which 
corresponds to a drop in the amplitude of the sixth 
harmonic by 15%. The amplitude of the first har¬ 
monic only drops by about 0.4%. Using a short 
segment reduces smearing of the pulse profile. (2) 
Typically 470 pulse profiles from the 300 s seg¬ 
ments are accumulated within each 4-day interval. 
Since we have a large number of profiles we can 
estimate sample variances for each Fourier coeffi¬ 
cient. Sample variances give us a measure of the 
noise present in the data and account for aperiodic 
noise if it is present. 

®In this paper, the words barycentered and barycentric refer 
to this correction. 
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The final step in our advanced pulsar search 
was a grid search in frequency using the set of 
typically several hundred estimated 20-50 keV 
pulse profiles from 4-day intervals of data. New 
XTE J1946-I-274 pulse frequencies were deter¬ 
mined from an initial grid search over 1200 evenly 
spaced trial barycentric frequencies in the range 
62.97-63.32 mHz for MJD 48361-51270 and 63.10- 
63.45 mHz for MJD 51270-51690. The search 
range was changed because continued spin-up of 
the pulsar moved the pulse frequencies outside 
the original search range. For each frequency off¬ 
set Aiy, a mean pulse profile is computed for each 
4-day interval. The mean pulse profile is of the 
form 

n 

m = X a/i exp(z27r/i(/)(t)) (1) 

h=l 

where 4>{t) = to), and 

tk is the midpoint barycentric time of segment k. 
The mean coefficients at are given by 

M 

^ Wkahk exp{-i2TrhA(j>k) (2) 

k=l 

with weights Wk = 

The typical statistic used for a search in pulse 
frequency or pulse period is the Z"^ statistic (Buc- 
cheri et al. 1983) given by 


= E 


\oth\ 


h=l 


(3) 


where crgh is the formal (i.e., Poisson statistical) 
error on the Fourier coefficient cth, respectively. 
If aperiodic noise is present, either due to XTF 
J1946-I-274 or due to other sources in the large 
BATSF field of view, e.g., Cygnus X-1, the Pois¬ 
son statistical error is an underestimate because 
the actual underlying noise level is higher than 
Poisson noise. This can create a problem with 
the statistic, causing it to depend on the noise 
level. To account for aperiodic noise, we used all 
of the typically ~ 470 pulse profiles from 300 s 
segments within each 4-day interval to estimate 
sample variances for the Fourier coefficients for 
each harmonic in the mean profile. Fach harmonic 
was treated separately. We then modified the 
statistic by replacing the Poisson variances with 
our computed sample variances to create a new 


statistic which we will call after Finger et al. 
(1999) where it is described in detail. Due to the 
large field-of-view of BATSF, other pulsars were 
often also present when we were measuring XTF 
J1946-I-274. If we limited our statistic to use the 
first 3 harmonics where XTF J1946-I-274 was the 
brightest, we reduced the chances of contamina¬ 
tion of the search results from other pulsars that 
happen to have harmonics near the higher har¬ 
monics of XTF J1946-I-274. The best-fit frequency 
for each 4-day interval was then determined using 
the Ya statistic. 

Using a similar method we also generated 
pulse frequency measurements for RXTE PCA 
observations during the initial outburst in 1998 
September-October and during the last 2 out¬ 
bursts, 2001 March-July. Barycentered Standard 
1 (125 ms, no energy resolution) data were fit¬ 
ted with a model consisting of a constant back¬ 
ground plus a 6th-order Fourier expansion in 
pulse phase model, creating an estimated 2- 
60 keV pulse profile for each PCA observation. 
The pulsed phase model consisted of a constant 
barycentric frequency, vq = 63.17466 mHz for 
the 1998 September-October observations and 
vq = 63.4239 mHz for the 2001 observations, esti¬ 
mated from projections of BATSF measurements. 
For each set of observations, we searched over a 
grid of 151 evenly spaced frequencies spanning the 
range r'oiO-17 mHz, the same size interval as used 
for BATSF. 

Figures la & c show the barycentered pulse fre¬ 
quency history and the 2-10 keV total flux his¬ 
tory, respectively, for XTF J1946-I-274. In Fig¬ 
ure la, the trend in the pulse frequencies repeat¬ 
edly changes from spin-up to spin down and back 
to spin-up at regular intervals (every 170 days). 
This repeated pattern is most likely due to the 
pulsar’s orbit. Within each 170 day cycle in the 
pulse frequencies, there are two outbursts. 

2.2. Orbit Fitting 

Examination of the observed pulse frequency 
history showed us two things: (1) a strong or¬ 
bital signature and (2) strong intrinsic torque vari¬ 
ations. To attempt to extract orbital parameters, 
we first fitted the observed pulse frequencies t'obs 
with a model consisting of a global orbit and a 
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global polynomial frequency model, i.e., 

^obs — ^emit{l /^) (4) 

where t'emit = Y^k=o^k{t — to)/kl is a model of 
the emitted frequency. The velocity relative to 
the observer P = v/c is given by 

27rax sint / (1 — cosujcosE — sinwsiniH 
PorbC \ {1-ecosE) 

_ (5) 

where ax sin i is the projected semi-major axis of 
the pulsar’s orbit; Porb is the orbital period; to is 
the periapse angle; and e is the orbital eccentricity; 

E is the eccentric anomaly, given by 

P — esinP = 27r'^-(6) 

and Tp is the epoch of periastron passage. We fit¬ 
ted models containing polynomials of orders 1-10. 
All of these models were very poor fits. Table 2 
lists the orbital parameters obtained using a first 
and 10th order polynomial for comparison. Errors 
on these parameters have been inflated by a fac¬ 
tor of (x^/d.o.f.)^/^, but are still expected to be 

considerably underestimated because the val¬ 

ues are so large. Figure 2 shows the emitted fre¬ 
quency model t'emit, ao(l ~ P) the constant term 
from the emitted frequency model times the ve¬ 
locity relative to the observer to illustrate the am¬ 
plitude of orbital effects, the full frequency model 
including the emitted frequency and orbit, and the 
frequency residuals for the 10th order polynomial 
model. 

Clearly a simple model such as a global poly¬ 
nomial did not suffice to describe the intrinsic 
torques. Instead we decided to model the intrinsic 
torques in a piecewise fashion. Each outburst was 
split into 2-6 intervals, depending on the length 
of the outburst, where each interval contained 
3-4 frequency measurements and spanned 16-24 
days or 16-32 days for the BATSE and the RXTE 
PCA data, respectively. In each interval, we fit¬ 
ted the observed pulse frequencies with a global 
orbit plus an independent linear frequency model. 
Equation 4 with t'emit = Vi + vpt — ti) as a model of 
the emitted frequency for each time interval i. Ta¬ 
ble 2 lists our best fit orbit. Figure 3 shows t'emit, 
J^i6(l — P) to illustrate the amplitude of orbital 
variations (Note was chosen because it lies in 
near the center of the frequency range.), the full 


frequency model including the emitted frequency 
and orbit, and the frequency fit residuals. This 
fit is not formally acceptable either; however, we 
believe this is mostly because this model still does 
not completely describe the intrinsic torques. Er¬ 
rors on the individual parameters in Table 2 have 
been inflated by (y^/d.o.f.)^/^. Comparing the 
three columns of Table 2 shows that the orbital 
parameters appear to be similar, whether we use 
a simple polynomial model or a more complicated 
piecewise model for t'emit- Figure lb shows the 
orbit-corrected spin frequencies. The regular pat¬ 
tern seen in Figure la has been removed, while 
intrinsic torque effects remain in the data. 

2.3. Orbital Phasing of Outbursts 

A close look at Figure 1 indicates that the out¬ 
bursts are not fixed in orbital phase and that we 
are seeing approximately 2 outbursts per orbit. To 
study the orbital phase of the outbursts, we fitted 
the 2-10 keV RXTE ASM flux in Figure Ic with 
a quartic polynomial, and used this polynomial to 
normalize the outbursts. We then computed the 
orbital position using 

X = a{cosE — e) (7) 

y = a(l - e^ji/^sinE (8) 

where a is the semi-major axis, estimated from 
ax sini assuming i = 70° which gives a companion 
mass of ~ 14.5M0. Figure 4 shows the normalized 
intensity versus orbital position, with darker colors 
indicating higher intensities. The outer loop of 
the spiral has a = 323 Rq, the estimated radius 
assuming i = 70°. For each successive orbit, we 
reduced a by 10% to offset the orbits. This plot 
spirals in solely to allow comparison of outburst 
phases. We are not claiming that the neutron star 
is spiraling in toward its companion. 

The outburst peaks clearly are not fixed in or¬ 
bital phase, although there is typically one peak 
on each side of the orbit. Just before periastron 
and just before apastron, the flux is typically rel¬ 
atively low. The cause of these low flux periods 
is unclear. Although the fluxes show considerable 
modulation, they remain detectable with RXTE 
throughout the period 1998 September - 2001 July. 
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2.4. Spectral Analysis 

To look for variations in energy spectra during 
the outbursts observed with the RXTE PCA, we 
first studied hardness ratios. Using FTOOLS® 5.1 
(Blackburn 1995) and Standard2 (16-s, 129 en¬ 
ergy channel) data, we created light curves in 4 en¬ 
ergy bands corresponding to 2-5, 5-10, 10-15, and 
15-20 keV for both the 1998 and 2001 outbursts. 
These light curves were background subtracted 
and their times were corrected to the Solar Sys¬ 
tem barycenter using FTOOLS. Next we formed 3 
hardness ratios from adjacent energy bands. Fig¬ 
ure 5 shows each of these ratios vs. total PCA 
source count rate in the 2-30 keV band. Figure 6 
shows color-color diagrams. In both figures, the 
grey-scale denotes intensity, with darker points de¬ 
noting higher intensities. (Note: This is the mea¬ 
sured 2-30 keV intensity, not the normalized 2-10 
keV intensity used in Section 2.3 and Figure 4.) 
For all of the hardness ratios, we see a correlation 
between hardness and intensity in the 2001 out¬ 
burst that is not present in the 1998 outburst. In 
the color-color diagrams, the hardness ratios are 
correlated for both sets of outbursts and have sim¬ 
ilar slopes. However, during the 2001 outbursts, 
points move from left to right along the correlation 
as intensity increases; while in the 1998 outburst, 
the points move around in hardness as intensity 
increases. 

To better quantify the observed spectral changes, 
we generated response matrices, count spectra, 
and background spectra using FTOOLS for PCA 
Standard2 data and HEXTE event mode data for 
each observation. We then fitted each of the ob¬ 
servations in XSPEC^ 11.1 (Arnaud 1996) with 
two models: (1) an absorbed power law with 
a high-energy cutoff and an iron line, XSPEC 
model: PHABS (POWERLAW -k GAUSSIAN) 
HIGHECUT^ and (2) an absorbed power law 
with an iron line and a cyclotron absorption fea¬ 
ture, XSPEC model: PHABS (POWERLAW -f 
GAUSSIAN) CYCLABS^. For both models we 
used only data in the 2.7-30 keV energy range for 
PCA data and in the 15-50 keV range for HEXTE. 
The power law normalizations were allowed to be 
independent for the PCA and each HEXTE clus¬ 
ter. Figure 7 shows spectral index, cutoff 

®http://heasarc.gsfc.nasa.gov/ftools 

^http://xspec.gsfc.nasa.gov 


energy, folding energy, and flux in the Fe line vs. 
2-60 keV flux. Figure 8 shows y^, spectral index, 
the cyclotron line energy, and the flux in the Fe 
line vs. 2-60 keV flux. Other spectral parameters 
did not show significant variations with intensity. 

A comparison of the two models showed that 
model 2 generally provided a better fit to the data 
(see Figures 7 & 8, top panel). In both mod¬ 
els, the power law photon index remains relatively 
independent of flux for fluxes above about 10“® 
ergs cm“^ s“^, however, below that flux level, the 
power law index softens as the intensity decreases. 
The integrated flux in the iron line is correlated 
with the 2-60 keV flux in both models, indicat¬ 
ing that the iron line is indeed a feature of XTE 
J1946-|-274’s spectrum and not a background fea¬ 
ture. In model 2 a cyclotron absorption feature, 
when it is constrained, tends to be near the value 
observed by Heindl et al. (2002). 

To estimate an upper limit to the quiescent flux 
from XTE J1946-I-274, we fitted PCA data from 
each of 6 observations after pulsations were unde¬ 
tectable in 2001 August - September with an ab¬ 
sorbed power law. These observations were consis¬ 
tent with a power law index of ~ 2.4, a hydrogen 
column of about 6x 10^^ cm“^, and a 2-30 keV flux 
of (4 — 7) X 10“^^ ergs cm“^ s“^. Since RXTE is 
not an imaging instrument, we cannot distinguish 
between low-level emission from XTE J1946-I-274 
and background sources such as the galactic ridge. 
However, this gives us an upper limit on the flux 
from XTE J1946-k274. 

2.5. Spin-up Torque-Flux correlations 

To estimate a bolometric correction to the 
RXTE ASM measurements, we plotted the 2-60 
keV fluxes from model 2 in Section 2.4 vs. the 
5-day average ASM count rate surrounding each 
PCA observation. We then fitted a line to these 
data, obtaining a bolometric correction of 1 ASM 
ct s“^ = (1.2 ± 0.1) X 10“® ergs cm“^ s“^(2-60 
keV). The error on this correction has been es¬ 
timated including systematic effects. Figure 9 
shows the best fit line to these data. Near the 
peak flux level, the ASM fluxes show evidence for 
a turnover relative to the PCA fluxes. We have 
attempted to account for this turn-over by includ¬ 
ing a systematic error of 10% on our bolometric 
correction. Using this bolometric correction, we 
then plotted the average bolometrically corrected 
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ASM flux vs. spin-up rate for each interval in 
our frequency model. Figure 10 clearly shows a 
correlation between spin-up rate and flux. Such a 
correlation suggests a disk is present because di¬ 
rect wind accretion is believed to be less efficient 
at transferring angular momentum (Ruffert 1997; 
Benensohn, Lamb, & Taam 1997; Ruffert 1999). 
If enough angular momentum is present, a disk 
will form. 

Simple accretion theory assumes that the ma¬ 
terial from the companion star is flowing onto a 
rotating neutron star with a strong magnetic field. 
This field determines the motion of material in a 
region of space called the magnetosphere. The size 
of this region is denoted by the magnetospheric 
radius r^, given by (Pringle & Rees 1972; Lamb, 
Pethick, & Pines 1973) 

Tm ^ (9) 

where G is the gravitational constant; Mx is 
the mass of the neutron star; and M is the 
mass accretion rate, which is assumed to be 
related to the observed bolometric flux F by 
M = 4ttcPFRx/GMx, where i?x is the radius 
of the neutron star and d is the distance to the 
pulsar. A: is a constant factor of order 1. Equa¬ 
tion 9 with k ~ 0.91 gives the Alfven radius ta 
for spherical accretion and with k ~ 0.47 n{ujs) 
gives the magnetospheric radius derived by Ghosh 
& Lamb (1979). 

The torque applied by accretion of matter onto 
a neutron star, assuming torques due to mat¬ 
ter leaving the system are negligible, is given by 
(Lamb, Pethick, & Pines 1973) 

^{2TTlv)=Mi ( 10 ) 


An accretion disk will form if the specific angu¬ 
lar momentum of the material accreted from the 
Be star’s disk is comparable to the Keplerian spe¬ 
cific angular momentum at the magnetospheric ra¬ 
dius, i.e., 

£ ~ C = {GMxrmf^^. (12) 

For XTE J1946-I-274, using the magnetic held 
measurement B = 3.1(1 + z) x 10^^ G (Heindl 
et al. 2002) which implies /r = 1.6 x 10^° G cm^, 
then £m = (2.7 — 2.8) x 10^^ cm^ s“^ for k = 0.91 
and Am = (1.9 — 2.0) x 10^^ cm^ s“^ for k = 0.47. 
For XTE J1946-f274, A = 0.9 - 1.3Am for k = 0.91 
and A = 1.3 — 1.9Am for k = 0.47, hence an ac¬ 
cretion disk is most likely present. In contrast, 
for the wind-fed system Vela X-1 where a disk is 
not expected to be present, z> ~ 6 x 10“^^ Hz 
s“^ (Inam & Baykal 2000), L ^ 2 x 10^® ergs 
s“^, and /i ~ 2.1 x 10®° G cm® (Makishima et 
al. 1999) leading to A ~ 3.5 x 10®"* cm^ s“® and 
Am cs 2.2 X 10®®" cm2 j ^ _ 0.002Am. Fur¬ 
ther, three-dimensional simulations of wind accre¬ 
tion show that the average specific angular mo¬ 
mentum accreted via wind accretion is always 
smaller than the Keplerian value (Ruffert 1997, 
1999). 

Since an accretion disk appears to be present, 
we can use simple accretion theory to derive a 
distance to XTE J1946-I-274. Substituting Equa¬ 
tions 12 into Equation 10 and solving for i> gives 

z> = 1.39fc®/2(GMx)”®/^Rx^^^”®Ai^/'^d®^/'^F°/'^. 

(13) 

Assuming k = 0.91 and the parameters defined 
earlier, 

z>i 2 = 5.8 X 10 - 242 ^^ 6 / 7 ^ 


where / is the moment of inertia of the neutron 
star and A is the specific angular momentum of 
the material. If I is assumed constant, then A is 
given by 

A = 27r/7'M-®, (11) 

where £> is the spin-up rate. To estimate A for 
XTE J1946-I-274, we assumed typical pulsar pa¬ 
rameters, Mx = IAMq, Rx = 10 km, I = 10®® g 
cm2, typical values of ~ 5 x 10-®2 Hz s”® and 
Fboi ~ 2 X 10“° ergs cm-2 s-®(see Figure 10), and 
a distance d = 8 — 10 kpc (Verrecchia et al. 2002). 
This yielded A = (2.4 — 3.8) x 10®®" cm2 g-i_ 


where h '12 is the spin-up rate in units of 10-®2 Hz 
s”® and Fg is the bolometric flux in units of 10“° 
ergs cm-2 g-i_ then fitted this model to the 
data in Figure 10. To account for the flux and i> 
errors, we computed x2 as 


X 


= E 


(F* - F“°‘®)2 - ;>niod('^mod))2 


(15) 

where F™”*® is a free parameter fitted to each 
flux and pmodj-^mod^ jg Equation 14 evaluated at 
F“°d. The solid line shows our best fit, with 
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X^/32 = 1.0, which gives us a distance of 9.5 ± 0.3 
kpc. This distance is consistent with optical obser¬ 
vations of the counterpart. The errors on the dis¬ 
tance are computed from the fit itself and take into 
account only the statistical errors on the flux and 
spin-up rate and the estimated systematic error on 
the bolometric correction. However, there is con¬ 
siderable uncertainty in our assumed neutron star 
parameters. The estimated distance depends on 
the following combination of neutron star param¬ 
eters: Using the neutron star 

equations of state given by Wiringa, Fiks, & Fab- 
riocini (1988) and assuming = 1.0 — I.SMq, 
we estimate the error on the distance due to un¬ 
certainty in the neutron star parameters is about 
30%, i.e., d = 9.5 ±2.9 kpc. 

2.6. Pulse Profiles 

Using FTOOLS and RXTE PC A event mode 
data, we created light curves and corresponding 
background light curves in 5 energy bands: 2-5, 
5-10, 10-15, 15-20, and 20-30 keV for each of 29 
observations, 12 from the 1998 outburst and 17 
from the 2001 outbursts. Each light curve was 
barycentered, background subtracted, and folded 
at the appropriate frequency from our frequency 
search. Each pulse prohle contained 50 phase bins. 
Eirst we examined prohles from the entire 2-30 
keV band to search for intensity dependent varia¬ 
tions. We aligned all of the pulse prohles by hnd- 
ing the phase of the minimum of the prohle using a 
quadratic interpolation and placing that minimum 
at phase 0.0. We found dehnite variations with in¬ 
tensity in both the 1998 and 2001 outbursts and 
also evidence that the pulse prohle was different 
at similar intensities during the rise and fall of the 
initial bright outburst. 

Prohles with similar shapes and similar intensi¬ 
ties were averaged to better illustrate the observed 
shape changes. Figure 11a shows the average 2-30 
keV pulse prohle at 6 different intensities. Fig¬ 
ure lib shows the peak-to-peak pulse fraction ver¬ 
sus mean hux for each prohle. In the following de¬ 
scriptions and in Figures 11a and b, the pulse pro¬ 
hles are numbered from one to six. (1) At the low¬ 
est intensities (6-17 counts s“^ PCU“^, 2-30 keV), 
during the 2001 outbursts, the prohle’s main fea¬ 
ture was a deep notch, which we used to align the 
prohles. The prohle consisted of an asymmetric 
structured main peak that peaked near phase 0.2. 


(2) As the intensity increased (23-65 counts s“^ 
PCU“^), a second notch becomes prominent near 
phase 0.45, making the prohle consist of two main 
peaks. These peaks were approximately equal in 
intensity at hrst, but the second peak brightens as 
the overall intensity increases. A small peak near 
phase 0.87 also appears. No RXTE PCA observa¬ 
tions were taken at intensities between 65 and 178 
counts s“^ PCU“^. (3) The next prohle shape oc¬ 
curs at intensities of 178-266 counts s“^ PCU“^, 
near the end of the 1998 outburst. This prohle 
is markedly different from that near the peak of 
the 2001 outburst (prohle 2). Both notches are 
broader, the peak near phase 0.65 is brighter than 
the peak near phase 0.3 and the small peak at 
phase 0.87 has disappeared. The main notch ap¬ 
pears to be less deep than at lower intensities. (4) 
Prohles at intensities of 285-308 and 330 counts 
s“^ PCU“^, all from the rise of the 1998 outburst, 
show a broader hrst peak and a less intense shoul¬ 
der following the second peak than the lower inten¬ 
sity prohles from the fall of the outburst (prohle 3). 

(5) At intensities of 322-328 counts s“^ PCU“^, 
during the decline of the 1998 outburst, the hrst 
peak is narrower than at similar intensities dur¬ 
ing the outburst rise (prohle 4), but broader than 
at lower intensities during the decline (prohle 3). 

(6) At the peak of the outburst, the width of the 
hrst peak is at an intermediate width between the 
prohles from the rise and fall of the outburst. Dur¬ 
ing the 2001 outbursts, the prohle shape appeared 
to depend primarily on intensity, while during the 
1998 bright outburst, the prohle shape depended 
on both intensity and whether the prohle was from 
the rise or fall of the outburst. 

Next we looked for variations in the pulse pro¬ 
hle versus energy. Figure 12 shows pulse prohles in 
5 energy bands from the peak of the 1998 outburst 
(left panel) and near one of the peaks in the pair of 
outbursts in 2001 (right panel). In the 1998 pro¬ 
hles, the prohle consists of two main peaks at lower 
energies (< 15 keV). The hrst peak is dominant at 
lower energies and the second peak becomes more 
dominant as energy increases. Also, as energy in¬ 
creases, the notch between the two main peaks hlls 
in. As the 2-30 keV intensity decreases, the pro¬ 
hle “shape” appears to move down in energy. For 
example, in the right hand panel of Figure 12, the 
2-5 keV prohle consists of two nearly equal peaks, 
reminiscent of the 10-15 keV prohle at the peak of 



the 1998 outburst. 

2.7. Optical Observations 

A regular monitoring program was established 
to study the strength and structure of the Ha opti¬ 
cal emission line during the X-ray observations in 
2001. This line arises from the circumstellar disk 
surrounding the Be star and the line properties are 
strongly related to the extent and dynamic struc¬ 
ture of the disk. Data were primarily collected 
from the 2.5m Isaac Newton Telescope (INT) in 
La Palma, Spain. The INT was equipped with the 
235-mm camera and EEV^IO CCD. The use of 
the R1200R grating results in a nominal dispersion 
of w 0.40A/pixel. Intermediate resolution spec¬ 
troscopy was also performed on 17th July 2001 and 
22nd October 2001 using the 1.93-m telescope at 
the Observatoire de Haute Provence, France. That 
telescope was equipped with the long-slit spectro¬ 
graph Carelec and the 1024 x 2048 EEV CCD. We 
used the 1200 In/mm grating in first order, result¬ 
ing in a nominal dispersion of « 0.45A/pixel. See 
Table 3 for details. 

The Ha profile obtained from each observation 
is presented in Figure 13. The lines have been 
shifted vertically by an arbitrary amount in order 
to present each profile clearly. The data of the 17 
July and 7 October 2001 are of a lower standard 
due to poor weather conditions. Nonetheless it 
is possible to accurately determine the equivalent 
width of each of the emission lines; values are pre¬ 
sented in Table 3. From these numbers one can see 
that there is little evidence for any changes in the 
line flux, assuming a constant continuum. How¬ 
ever, from Figure 13, one can see obvious shifts in 
the line position after 2001 July. 

If one takes the first spectrum (3 May 2001) 
as a base line and subtracts it from each of the 
other spectra after they have been normalized to 
the same peak value, then significant changes are 
apparent. The result of this process is shown in 
Figure 14. It is apparent that little happens to the 
profiles (and, presumably, the circumstellar disk) 
until July/August 2001. Starting with the spec¬ 
trum of 17 July 2001 one can see significant per¬ 
turbations occurring in the line profile indicative 
of density changes in the disk structure. Between 
August and September the direction of the per¬ 
turbation changes, suggesting that the density en- 
hancement/rarefication had changed sides (or ro¬ 


tated around) in the disk. The size of the per¬ 
turbation also seems to have been increasing from 
July to October, though the separation of the red 
and blue peaks remains approximately constant at 
~340 km/s. 

3. Discussion 

XTE J1946-I-274 is a very unusual Be/X-ray 
binary system and its behavior does not fit well 
into the standard normal/giant outburst behav¬ 
ior. Instead it showed an extended period of ac¬ 
tivity from 1998 September - 2001 July when the 
X-ray flux showed considerable modulation, with 
two peaks or outbursts per orbit, but the X-ray 
flux never dropped below the detection threshold 
of RXTE. This extended period of activity more 
closely resembles a series of normal outbursts than 
a single giant outburst; however, unlike typical 
normal outbursts, these outbursts shift rapidly in 
orbital phase, there are two outbursts per orbit, 
and the X-ray flux does not drop dramatically 
between outbursts. The initial outburst shows 
considerable spin-up v > 10“^^ Hz s“^, like a 
giant outburst, but this spin-up is not substan¬ 
tially larger than that seen in later outbursts. In 
addition, the peak flux (2-60 keV) of the initial 
outburst is about 7 x 10“® ergs cm“^ s“^, while 
the peak flux (2-60 keV) of the last two outbursts 
observed with the PCA is about 2 x 10“® ergs 
cm“^ s“^, only a factor of ~ 3.5 fainter than the 
brightest outburst. These fluxes correspond to lu¬ 
minosities of 8 X 10^^ ergs s“^and 2 x 10^^ ergs 
s“^, respectively, assuming our best-fit distance of 
9.5 kpc. 

In giant outbursts of Be/X-ray binaries, ac¬ 
cretion disks are expected to be present and in¬ 
deed, evidence for an accretion disk, based on cor¬ 
relations between the observed flux and spin-up 
rate, has been found for several sources during gi¬ 
ant outbursts (Parmar et al. 1989; Reynolds et 
al. 1996; Bildsten et al. 1997; Wilson et al. 1998; 
Stollberg et al. 1999). Independent evidence for 
an accretion disk based on the detection of quasi- 
periodic oscillations during a giant outburst has 
been found for EXO 2030-1-375 (Angelini, Stella, 
& Parmar 1989) and A0535-I-262 (Finger, Wilson, 
& Harmon 1996). Until recently, normal outbursts 
were believed to be due to direct wind accretion 
from the Be disk, so significant spin-up was not 
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Equation 9 and the corotation radius is given by 
rco = (GM)1/3(2^:.)-2/3 ( 16 ) 

where v is the spin frequency of the pulsar. Setting 
T-m = ?'co gives the threshold flux for the onset of 
centrifugal inhibition of accretion, i.e., 


expected because direct wind accretion is not be¬ 
lieved to be very efficient at transferring angular 
momentum (Ruffert 1997; Benensohn, Lamb, & 
Taam 1997; Ruffert 1999). If enough angular mo¬ 
mentum is present in the accreted material, an ac¬ 
cretion disk will form. However, evidence for spin- 
up during normal outbursts has been observed in 
GS 0834-430 (Wilson et al. 1997), 2S 1417-624 
(Finger, Wilson, & Chakrabarty 1996), 2S 1845- 
024 (Finger et al. 1999), and in EXO 2030-f375 
(Stollberg et al. 1999; Wilson et al. 2002). 

We see a correlation between spin-up and flux 
for XTE J1946-I-274 as shown in Figure 10. The 
spin-up rate and its correlation with bolometric 
flux during XTE J1946-|-274’s outbursts suggest 
an accretion disk may be present. In addition, 
our calculations show that the specific angular mo¬ 
mentum of the accreted material is comparable to 
the Keplerian specific angular momentum, hence a 
disk is expected to form in XTE J1946-I-274. Fur¬ 
ther, our fit to the spin-up vs. bolometric flux 
correlation, which assumed a disk was present, 
yielded a distance of 9.5 ± 2.9 kpc which is consis¬ 
tent with the distance of 8-10 kpc derived from op¬ 
tical observations Verrecchia et al. (2002). Hence 
an accretion disk is likely to be present in XTE 
J1946-f274. 

Between 2001 July 31 (MJD 52121) and 2001 
August 9 (MJD 52130), XTE J1946-f274 dropped 
below the PCA’s detection threshold. At approx¬ 
imately the same time, between 2001 June 29 
(MJD 52089) and 2001 July 17 (MJD 52107), the 
Ha profile began to change rapidly. The profile 
was stable for 3 observations on 2001 May 3, May 
10, and June 29 (MJD 52032, 52039, & 52089) that 
corresponded to X-ray observations during the de¬ 
cline of the second to last outburst, during the low 
state between the last two outbursts, and near the 
peak of the last outburst. The coincidence of the 
change in Ha with the X-ray turn-off suggests that 
changes in the Be disk caused the X-ray outbursts 
to cease. 

To determine whether or not the observed X- 
ray turn-off was due to centrifugal inhibition of 
accretion (Stella, White, & Rosner 1986), we esti¬ 
mate the flux at the onset of this effect by equat¬ 
ing the magnetospheric radius to the corotation 
radius. The magnetospheric radius is given by 


cs 3.0x10 ® ergs cm ^ g ^ 2 / 3 ^^ 

(17) 

where /igo, Mi. 4 , Rg, and Pis.ss are the pulsar’s 
magnetic moment in units of 10^° G cm^, mass in 
units of 1.4 M 0 , radius in units of 10® cm, and 
spin period in units of 15.8 seconds, respectively. 
Using d = 9.5 kpc and = 1.6 x 10®® G cm® 
yields ci 1.2 x 10“®® ergs cm“^ . For 

fc = 0.47 - 0.91, E™ ~ (0.6 - 6.0) x 10“®® ergs 
cm“^ s“®. Our measured upper limit fluxes are in 
the range (4 — 7) x 10“®^ ergs cm“^ s“®, consis¬ 
tent with XTE J1946-I-274 entering the centrifugal 
inhibition of accretion regime when the X-rays be¬ 
came undetectable with the PGA. 


We see two outbursts per orbit in XTE J1946-I-274; 
however, these outbursts are not fixed in orbital 
phase. This combination of behaviors is unique to 
XTE J1946-I-274 and is likely not due to a single 
mechanism. The optical data and mass function 
suggest that the orbital plane and the Be disk are 
not aligned in this system. Such a misaligned sys¬ 
tem would be expected to produce two outbursts 
per orbit, each approximately corresponding to 


4. Conclusions 

We propose that XTE J1946-I-274 is a system 
in which the Be star’s equatorial plane and the 
orbital plane are not aligned. This is suggested 
by: (1) The optical observations show a relatively 
narrow FWHM (8.6 A on 2001 May 10, which cor¬ 
responds to ± 200 km s“®) for the single-peaked 
Ha line, indicating that the Be star is viewed from 
a relatively low inclination angle, i.e., nearly pole- 
on. (2) The orbital signature is quite obvious in 
the pulse frequencies, indicating that the orbital 
inclination angle is likely not low: the derived 
mass function f{M) = from the piece- 

wise frequency model indicates an inclination an¬ 
gle > 46° for the expected mass range of 10—I 6 M 0 
derived from optical observations. Figure 15 shows 
the range of allowed inclinations using the 1 -ct er¬ 
rors on the mass function. 
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the neutron star’s passage through the Be disk; 
hence these outbursts would be expected to be 
fixed in orbital phase. For XTE J1946+274, we 
are faced with the difficult task of explaining not 
only why we see the expected two outbursts per 
orbit, but also why these outbursts are not fixed 
in orbital phase. We propose that given the rel¬ 
atively high luminosities and spin-up rates in all 
of the outbursts that perhaps XTE J1946-|-274’s 
unusual outburst behavior can be explained as a 
giant outburst. In normal outbursts, the Be disk 
is expected to be truncated at a resonance ra¬ 
dius by tidal forces from the neutron star’s orbit 
(Okazaki & Negueruela 2001); however, in giant 
outbursts, where much more material is believed 
to be present in the disk, it is believed that the 
disk is no longer truncated. Hence more mate¬ 
rial likely means a larger disk around the Be star. 
Density perturbations propagating in the Be disk 
could change the contact points of the disk and 
the orbit. These perturbations would have to be 
moving quite rapidly in the disk. This idea of a 
giant outburst combined with density perturba¬ 
tions seems most likely to work in a system where 
the inclination angle between the Be disk and the 
orbital plane is fairly small. It is not clear that 
this angle is small in XTE J1946-I-274. We do not 
claim that the idea of a giant outburst in an in¬ 
clined system fully explains the behavior observed 
from XTE J1946-I-274. Instead we put the idea 
forth as a suggestion to those doing simulations of 
these systems. 

The phasing of the outbursts of XTE J1946-I-274 
(see Figure 4) is difficult to understand. How¬ 
ever, it is not the first Be/X-ray binary to show 
shifts in outburst phase. EXO 2030-1-375 under¬ 
went a fairly sudden (within 4 orbits) shift in 
the peak phase of the outbursts from 6 days af¬ 
ter periastron (phase ~ 0.13) to 2.5 days before 
periastron (phase -0.05), followed by a gradual 
shift in outburst phase to 2.5 days after perias¬ 
tron (phase 0.05). This shift was believed to be 
associated with a density perturbation observed 
in the Be disk via Ha observations (Wilson et 
al. 2002). GS 0834-430 also underwent sudden 
shifts in outburst phase, undergoing 9 outbursts 
approximately centered on periastron, followed by 
dramatic shifts to outbursts centered on phases 
0.37 and 0.75 (Wilson et al. 1997). A companion 
was not known at the time for GS 0834-430, so 


we cannot confirm that density perturbations in 
the Be disk were responsible for these phase shifts. 
In XTE J1946-I-274, we have a similar problem, 
with the companion not observed for much of the 
time the X-ray outbursts were occurring. In addi¬ 
tion, the problem is compounded by the fact that 
we are likely viewing the Be star nearly pole-on; 
hence the range of projected rotational velocities is 
very small, making it difficult to detect perturbed 
velocities which indicate density perturbations. 
Only large density perturbations, such as the one 
that coincided with the X-ray turn-off, are likely 
to be detected. 

We are grateful to Luisa Morales for provid¬ 
ing one of the optical spectra. This research has 
made use of data obtained from the High En¬ 
ergy Astrophysics Science Archive Research Gen- 
ter (HEASARG), provided by NASA’s Goddard 
Space Flight Center (GSFC). RXTE ASM quick- 
look results were provided by the ASM/RXTE 
teams at MIT and at the GSFC SOF and GOF. 
The INT is operated on the island of La Palma by 
the Isaac Newton Group in the Spanish Observa- 
torio del Roque de los Muchachos of the Institute 
de Astrofisica de Canarias. Based in part on obser¬ 
vations made at Observatoire de Haute Provence 
(CNRS), France. 
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Fig. 1.— XTE J1946-I-274 outburst history, (a) 
Barycentered pulse frequency measurements with 
BATSE (filled circles) and RXTE PCA (open 
circles), (b) The pulse frequency measurements 
shown in (a) have been corrected for our orbital 
solution listed in the fourth column of Table 2 (c) 
2-10 keV flux history measured with the RXTE 
ASM. Each point is a 10-day average flux. Filled 
squares denote 3-a or better detections. The solid 
curve denotes a quartic polynomial fitted to the 
RXTE ASM flux (see Section 2.3). Vertical dot¬ 
ted lines in all three panels denote the estimated 
times of periastron passage using the orbital pa¬ 
rameters in the fourth column of Table 2. 


This 2-column preprint was prepared with the AAS lAT^ 
macros v5.0. 
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Table 1 

Observations of the Sky Region Including XTE J1946+274 


Dates 

MJD 

Instrument 

No. Obs.'‘ 

Total PCA Exposure (ks)^ 

1991 Apr 15 - 2000 May 27 

48361-51691 

BATSE 



1996 Feb 23 - 2002 Apr 11 

50136-52376 

RXTE ASM 



1998 Sep 16 - 1998 Oct 14 

51072-51100 

RXTE PCA, HEXTE 

12 

85.87 

2001 Mar 9 - 2001 Sep 25 

51977-52177 

RXTE PCA, HEXTE 

24 

129.49 


^For pointed observations. 


Table 2 
Orbit Fit 


Parameter 

Linear i/emit 

10th Order i^emit 

Piecewise i/emit 

Porb (days) 

173 ± 2 

167.8 ± 0.6 

169.2 ± 0.9 


JD2451524.6 ± 16.9 

JD2451572.3 ± 1.6 

JD2451558.7 ± 4 

ax sin i (It-sec) 

972 ± 118 

474 ± 14 

640 ± 120 

e 

0.19 ± 0.13 

0.35 ± 0.03 

0.33 ± 0.05 

<jj (degrees) 

-153 ± 37 

-38 ± 5 

-91 ± 23 

/(M) (Mq) 

32 ± 11 

4.0 ± 0.3 

Q y+6.9 
' -4.5 

to 

p. 

0 

2699556/101 

28756/92 

220.1/37 


Table 3 

Journal of optical observations. 


Date 

Telescope^ 

Instrument^ 

Grating 

Detector 

Ha EW (A) 

3 May 2001 

INT 

IDS+500 camera 

R1200R 

TEK5 

45.9±0.6 

10 May 2001 

INT 

IDS+235 camera 

R1200R 

EEVIO 

41.0±0.8 

29 Jun 2001 

INT 

IDS+235 camera 

R1200R 

EEVIO 

41.1±0.3 

17 Jul 2001 

OHP 

Carelec 

12001/mm 

EEV 

39.1±0.4 

12 Aug 2001 

INT 

IDS+235 camera 

R1200R 

EEVIO 

43.9±0.3 

26 Sep 2001 

INT 

IDS+235 camera 

R1200R 

EEVIO 

43.6±0.3 

7 Oct 2001 

INT 

IDS+235 camera 

R1200R 

EEVIO 

38.5±0.5 

22 Oct 2001 

OHP 

Carelec 

12001/mm 

EEV 

42.5±0.3 


^INT = Isaac Newton Telescope, OHP = Observatoire de Haute Provence 
^IDS = Intermediate Dispersion Spectrograph. 
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Fig. 2.— Top: The solid curve denotes the emit¬ 
ted frequency modeled as a 10th order polyno¬ 
mial. The dotted curve illustrates the magnitude 
of orbital effects. The open circles denote the full 
pulse frequency model including both the 10th or¬ 
der polynomial and the orbital effects. Bottom: 
Pulse frequency residuals for the model consisting 
of a global orbit and a global 10th order emitted 
frequency model. 



Li_ 


51000 51400 51800 52200 

Time (MJD) 

Fig. 3.— Top: The broken curve denotes the 
piecewise emitted frequency model. The dotted 
curve illustrates the magnitude of orbital effects. 
The open circles denote the full pulse frequency 
model including both the piecewise emitted fre¬ 
quency model and the orbital effects. Bottom: 
Pulse frequency residuals for the model consisting 
of a global orbit and piecewise emitted frequency 
model. 





Periastron 


Fig. 4.— Normalized intensity (darker colors in¬ 
dicate higher intensities) vs. position in the orbit. 
The outside loop is the shape of the actual or¬ 
bit. The semi-major axis of each successive orbit 
is reduced by 10% to offset the orbits, allowing 
comparison of the orbital phase of the outburst in 
successive orbits. This plot is not meant to indi¬ 
cate that the neutron star is spiraling in toward 
its companion. 
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Fig. 5.— Hardness-intensity diagrams. Filled 
squares denote points from the 1998 outburst and 
filled circles denote the 2001 outbursts. The color 
darkens as the intensity increases within each out¬ 
burst to allow comparison with Figure 6 
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Fig. 6.— Color-color diagrams. Filled squares de¬ 
note points from the 1998 outburst and filled cir¬ 
cles denote the 2001 outbursts. For both sets of 
outbursts, the color darkens as the intensity in¬ 
creases. 



Fig. 7.— Spectral fitting results for an absorbed 
power law with a high energy cutoff and a Gaus¬ 
sian iron line. Shown are the goodness of fit 
X^/d.o.f., the power law photon index, the cut off 
energy Flcut, and the integrated flux in the iron 
line versus 2-60 keV flux. Filled squares denote 
the 1998 outburst and open squares denote the 
last pair of outbursts observed in 2001. 



Fig. 8.— Spectral fitting results for an absorbed 
power law with a cyclotron absorption feature and 
a Gaussian iron line. Shown are the goodness of 
fit x^/d.o.f., the power law photon index, the cy¬ 
clotron line energy Ecyc, and the integrated flux 
in the iron line versus 2-60 keV flux. The dotted 
line denotes the cyclotron energy measured during 
the 1998 outburst by Heindl et al. (2002). Filled 
squares denote the 1998 outburst and open squares 
denote the last pair of outbursts observed in 2001. 



Fig. 9.— ASM count rate vs. PCA 2-60 keV flux 
measurements. The solid line denotes our best fit 
the the correlation. This fit was used to estimate 
a bolometric correction to the ASM data. Dashed 
lines correspond to our best fit slope ±10% in¬ 
dicating our estimate of systematic error due to 
changing beaming angles at high fluxes. 
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Fig. 10.— Spin-up rate vs. bolometrically cor¬ 
rected ASM flux. The solid line is our best fit of 
the simple torque model with a distance of 9.5 kpc. 
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Fig. 11.— (a): Average 2-30 keV pulse profiles 
from 6 intensity bands versus intensity. Profiles 
are numbered 1-6, corresponding to descriptions 
in the main text. Profiles 1 and 2 are from the 
2001 outbursts. Profiles 4 and 5 are from similar 
intensities during the rise and fall, respectively, of 
the 1998 outburst. Each profile has been normal¬ 
ized to its mean flux. Profiles are offset by ar¬ 
bitrary amounts to keep them from overlapping, 
(b): Peak-to-peak pulse fraction versus mean 2-30 
keV flux for the pulse profiles shown in panel (a). 
Horizontal bars indicate the flux range spanned by 
the profiles included in the average. 
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Fig. 12.— Pulse profiles in 5 energy bands (from 
top to bottom) 2-5, 5-10, 10-15, 15-20, 20-30 keV 
from the peak of the 1998 outburst (left) and the 
2001 outbursts (right). 


H al])lia s])ectra of XTF J1946 + y7'l 



Fig. 13.— Ha profiles from 7 observations of XTE 
J1946-F274 during our 2001 X-ray observations. 
The small dip near 6470 A in the 10 May 2001 
spectrum is most likely due to a cosmic ray hit. 
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Fig. 14.— Ha profiles in Figure 13 were normal¬ 
ized to the same peak value and then the profile 
from 2001 May 3 was subtracted from all other 
profiles to illustrate differences in the profiles. At 
some time between 2001 Jun 29 and Jul 17, the 
profile began to change rapidly. 



COS i 


Fig. 15.— Companion Mass versus the cosine of 
the inclination angle. The solid line is the mass 
function from the fourth column of Table 2. Dot¬ 
ted lines denote the 68% confidence error region 
on the mass function. Dashed lines indicate the 
range of companion masses derived from optical 
observations. The hatched region denotes the al¬ 
lowed masses and inclinations. 
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